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EFFECTS OF THERMAL LOADING ON COMPOSITES WITH CONSTITUENTS 
OF DIFFERING THERMAL EXPANSION 
by Charles A. Hoffman 
Lewis Research Center 

SUMMARY 


A heuristic study was conducted to analytically and experimentally evaluate the pos- 
sible effects of differences in coefficients of thermal expansion of composite constituents 
on the structural integrity and performance of metal-metal composites. Both laminated 
sheet or foil and fiber types of composites were considered, and special reference was 
made to refractory- metal- reinforced superalloy matrix composites. Approximate equa- 
tions were derived to permit the calculation of elastic stresses and elastic-plastic 
stresses and strains. 

Stresses were calculated for model tungsten- reinforced - 80Ni+20Cr matrix com- 
posites of both laminate and fiber configuration. Heating and cooling between 80° and 
2000° F (26.5° and 1093.5° C) was assumed. The elastic tensile, compressive, and 
peak shear stresses on heating or cooling could be far in excess of the tensile strengths 
and estimated shear strength. Approximate elastic -plastic analysis resulted in lowered 
stresses, but indicated concurrent strains of sufficient magnitude to cause thermal fa- 
tigue if heating and cooling cycles were encountered. 

Limited experimental studies with tungsten-reinforced - 80Ni+20Cr- matrix laminate 
composites and slow heat-cool cycles between room temperature and 2000° F (1093. 5° C) 
caused observable structural damage in 3 and 11 cycles for laminates with 0.020 and 
0.005 inch (0.05 and 0.0125 cm) lamina, respectively; a 0.001-inch (0.0025-cm) 
lamina composite appeared undamaged at 11 cycles. Rapid heat-cool cycles produced 
observable damage in one, three, and six cycles for composites with 0.020-, 0.005-, 
and 0. 001-inch (0. 050- , 0.0125-, and 0. 0250-cm) lamina, respectively. Six heat- cool 
cycles between room temperature and 1600° F (877° C) caused external change and 
microfracture of the matrix in a 0. 70 volume fraction tungsten-fiber - copper-matrix 
composite; a 0.40 volume fraction tungsten-fiber composite sustained an external ma- 
trix crack after six such cycles. Heating the 0. 70 volume fraction tungsten-fiber com- 
posite from room temperature to 1652° F (900° C) in a hot-stage microscope produced 
microfracture in the matrix. A 0.40 volume fraction tungsten-fiber composite exhibited 
severe matrix strain as a result of being heated from room temperature to 1652° F 
(900° C). 



INTRODUCTION 


There is, currently, considerable interest in composite materials: such materials 
offer many potential and demonstrated advantages for achieving higher use-temperatures 
and strength to weight ratios. Composites with superalloy matrices and refractory fiber 
reinforcements are being extensively studied in an effort to obtain materials for use at 
elevated temperatures (refs. 1 to 3). Some combinations of materials offer substantial 
promise for use in creep-rupture applications in the 2000° to 2200° F (1093. 5° to 
1205° C) range (ref. 2). Numerous other metal-fiber matrix composites, such as boron 
in aluminum, boron in titanium, boron in magnesium, silicon carbide in aluminum, sili- 
con carbide in magnesium, and silicon carbide in titanium are also being considered, 
but for lower use-temperatures. Laminated sheet or foil composites have also received 
attention as a means of achieving more impact resistant or fracture resistant materials 
(refs. 4 to 6). Laminate composites can comprise alternating strong reinforcement 
laminae and ductile matrix laminae, bonded together. Laminated composites presum- 
ably can be relatively simply fabricated and may have less anisotropy than fiber- 
reinforced composites. 

The successful use of dissimilar materials when used in composites, requires ade- 
quate compatibility between the materials. The complimentary aspects of a materials 
combination must outweigh any deleterious aspects. Many combinations of materials 
considered for use in reinforced metal-matrix composites have appreciable differences 

in the coefficients of thermal expansion of the reinforcement and the matrix. Typical 

— 

superalloy matrices have coefficients of thermal expansion of about 9x10 to 10x10" 
inch per inch per °F (16x10”® to 19x10“® m/(m)(°C)) over the temperature range of 80° 
to 2000° F (26.5° to 1093. 5° C). The refractory metals that might be used for rein- 
forcement (and presumably their alloys also) have coefficients of thermal expansion of 
about 3x10 D inch per inch per °F (54x10” m/(m)(°C)), over the same temperature 
range. The aluminum, titanium, and magnesium matrices can also have appreciable 
differences in their coefficients of thermal expansion compared with those of their rein- 
forcements. 

One would expect that thermal treatment (i. e. , temperature cycling) given to mate- 
rials having sufficient differences in thermal expansions could damage the composite by 
disrupting its structural integrity. Adverse effects could occur on cooling (or heating) 
the composite after initial consolidation treatment or during use when heating and cool- 
ing cycles are part of the exposure conditions. 

This possible problem that could arise from the generation of residual thermal 
stresses has been recognized by others (e. g. , refs. 3 and 7 to 9). Scala (ref. 7) calcu- 
lated that elastic tensile stresses (i. e. , 100 to 200 ksi (689 to 1380 MN/m^)) far above 
matrix yield and ultimate tensile strengths could be developed in the matrices of some 
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composites. He also suggested that composites might fail as a result of shear stress 
along the bond interfaces. In reference 9, an analysis was made to determine the limit- 
ing coefficients of thermal expansion of possible metal reinforcing fibers incorporated 
into a ceramic matrix (e. g. , ZrOg). These limiting coefficients would then prevent ma- 
trix damage from occurring as a result of differential expansions due to assumed tem- 
perature changes. 

The purpose of the present study was to gain a further understanding of the possible 
effects of differences in thermal expansions of metal-metal composite constituents on 
composite integrity and performance. This heuristic study has the following specific 
obj ectives : 

(1) To develop approximate mathematical expressions for the internal stresses that 
might arise and to determine the approximate magnitudes of these stresses and their as- 
sociated strains for a representative foil laminate composite and a representative fiber 
composite. 

(2) To obtain experimental indication of the effect of thermal cycling and resultant 
thermal stress on composite integrity. 

(3) To use the concepts gained from the analytical results and the experimental study 
results and to consider the practical significance of internal stresses and strains and 
ways of meliorating any possible adverse effects. 

The scope of this study includes the derivation of approximate equations for elastic 
stresses that might be generated in laminate composites and in fiber composites because 
of thermal stresses arising from differences in coefficients of thermal expansion. The 
equations and auxiliary methods also permit the estimation of elastic-plastic stresses 
and strains. In addition, limited experimental studies are made on laminate and fiber 
composites to obtain practical indications of the possibility of composite damage from 
thermal stress-strain effects due to thermal expansion differences. The mathematical 
analysis is not intended to be rigorous at all times; such a rigorous treatment is beyond 
the scope of this report. Where existing derivations or methods of analysis were appli- 
cable, they were used in this study. 


SYMBOLS 

A total area 
a area of fiber of sheet 

b fiber radius 

c radius of matrix associated with fiber (see appendix D) 

D diameter of fiber (i. e. , reinforcement) or associated matrix 
E modulus of elasticity 
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G shear modulus 

K ratio of transverse to longitudinal stresses (see appendix E) 
l length of sheet or fiber 

m matrix breadth 

N number of elements (fibers or laminae) 

P total force or load 

p pressure 

r reinforcement breadth 

T temperature 

t thickness 

u matrix displacement 

V volume fraction 

W width of sheet 

x distance in x direction (see subscripts) 

a coefficient of linear thermal expansion 

0 (b/c) 2 

y shear angle 

e mechanical strain 

X parameter, see appendixes B and E 

v Poisson's ratio 

P (r/c) 2 

cr tensile or compressive stress 

r shear stress 

Sul (scripts: 
av average 

c composite 

E elastic 

e equivalent 

h homogeneous 

1 internal 

m matrix 

o external 
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plastic 

radial direction 
reinforcement 
total mechanical strain 
orthogonal directions 
tangential direction 

subscripts differentiating a given symbol 

LAMI NATE COMPOSITES 

STRESSES IN LAMINATE COMPOSITES MADE FROM SHEET OR FOIL 

Assumptions 

The following assumptions are part of the analysis : 

(1) Constituents are perfectly bonded to each other. 

(2) Each constituent is homogeneous and isotropic, and each has the same Poisson's 
ratio. 

(3) Tensile or compressive stresses are uniform through the thickness of each con- 
stituent. 

(4) The loading is such that a = a: that is, a biaxial stress state exists in the 

x y 

laminae. 

(5) Laminae thickness is assumed to be small compared with the other dimensions; 
therefore, <j z is considered to be equal to zero. 

(6) No interphase region exists. 

Edge effects are not considered. 


P 

R 

r 

t 

x, y,z 

e 

1,2 


Mechanism of Stress Generation 

The model is illustrated in figure 1. It is composed of two different sheet or foil 
materials, each with a different coefficient of thermal expansion. These two materials 
are alternately positioned and bonded together. It is assumed for this analysis that 
stress-free bonding has been accomplished at elevated temperature and that the compos- 
ite is subsequently cooled down to room temperature. The same equation would be de- 
rived if stress-free bonding could be accomplished at room temperature and the compos- 





Figure 2. - Forces acting on element of laminate composite 
at some point away from ends; after cooling. 


Figure 1. - Geometry of laminate composite showing alternating matrix 
and reinforcement laminae and element considered in the shear analysis. 


ite heated. It is apparent that, because of the bond between constituents, the constituent 
with the higher coefficient of thermal expansion will sustain a residual tensile stress (in 
the x and y directions) on cooling, and the constituent with the lower coefficient of ther- 
mal expansion will sustain a compressive stress (in the x and y directions) on cooling 
as indicated in figure 2. The reverse would be true if bonding took place, without stress, 
at room temperature and if the composite were then heated to an elevated temperature. 
The biaxial stress state (assumption (4)) will produce an equivalent stress. (Subsequent 
heating or cooling obviously would also produce stresses and strains, but not necessarily 
the same as those which occur on the initial thermal cycle; this because of possible 
changes in material properties or stress state. ) The actual derivations are given in ap- 
pendix A. The resultant equations for equivalent stress are 


V E E 
rrm 


em 


(V E + V E )(1 - v) L 
' r r m m A ’ 


(« r “ “m )AT + - teonn - J 


erp emp' 


( 1 ) 


°er 


V E E 
m m r 


(V E + V_E )(1 - v ) 
v » r ^ r T * m m /v ’ 


(« m “ a r )AT (e 


emp erp 


( 2 ) 
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These equations may be used to calculate elastic and elastic-plastic stresses as well as 
elastic-plastic strains. 


ELASTIC SHEAR STRESSES IN LAMINATE COMPOSITES 

Assumptions 

(1) Reinforcement and matrix are perfectly bonded to each other and are homoge- 
neous, isotropic, and have the same Poisson's ratio. 

(2) A state of equal biaxial stress exists (i. e. , ct x = a ). 

(3) The stresses are in the elastic range. 

(4) The stress a x is independent of y, away from edges. 

(5) The thickness of the matrix is equal to or less than that of the reinforcement. 

(6) No interphases exist. 

Interacting edge effects are not considered (i. e. , an element as shown in fig. 1 is con- 
sidered). 


Mechanism of Stress Generation 

As a result of either heating or cooling a laminate composite after stress-free low- 
or high-temperature consolidation, shear stresses will be developed in the matrix of the 
composite since the matrix will expand or contract more than the reinforcement. The 
stresses developed by shear in the matrix after cooling and the strains developed as a 
result of matrix shear after cooling are illustrated in figures 3 and 4, respectively. 
These figures are discussed more fully in the derivation of the shear-stress equation in 
appendix B. The resultant equation for elastic shear stress is 


where 


T mzx 



(a -a )AT 
v m r J 

» - y m> / E r\ 

< X - "r> VV 



-»A m 

2G e 
m 


X 



“'W 1 - 


E t 
r r 


(3) 


( 4 ) 
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X 

(a) Typical element (b) Forces on reinforcement element. 

Figure 3. - Element of laminate composite being considered and forces 
acting on the reinforcement after cooling. 


X 



Line of symmetry 
within matrix 

Figure 4. - One-half thickness of matrix element (side view, near tip) 
before cooling (initial position) and after cooling (final position). 
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CALCULATION OF RESIDUAL THERMAL STRESSES AND STRAINS IN 
MODEL SHEET OR FOIL LAMINATE COMPOSITES 


Residual thermal stresses and strains have been calculated for sheet or foil rein- 
forced laminate composites using equations (1) to (3). The calculations were made for 
two model specimens: one was assumed to have been consolidated at 2000° F (1093. 5° C) 
(stress free) and then cooled to room temperature (condition A); the other consolidated 
in a stress-free manner at room temperature and heated to 2000° F (1093. 5° C) (condi- 
tion B). Condition A would also be analogous to a specimen exposed to 2000° F 
(1093.5° C) during use and cooled to room temperature; condition B would be analogous 
to a specimen being at room temperature and then exposed to 2000° F (1093. 5° C). In 
both instances, the laminate composite specimens were assumed to be stress free at 
their respective initial temperatures. In these calculations the reinforcement is tung- 
sten and the matrix is Nichrome V (80Ni+20Cr). The corresponding properties of im- 
portance are presented in table I. 

The purpose of the calculations presented herein is to obtain indications of the mag- 
nitudes of thermally induced stresses and strains that might occur in use for tungsten 
reinforced 80Ni+20Cr. The calculated values of stresses and strains would be compara- 
ble to those in many other potential refractory-metal-reinforced - superalloy- matrix 
composites since the respective coefficients of linear thermal expansion are similar as 
are the moduli, Poisson's ratios, and yield strengths (see tables II to IV). The elevated 
temperature elastic modulus for the 80Ni+20Cr may be lower than that for other possible 
matrix materials so that those other materials could generate larger thermal stresses 
on heating. 

TABLt L - PROPERTIES OF CONSTITUENTS IN MODEL LAMINATE COMPOSITE 


(Poisson's ratio is assumed to be 0. 25 for both materials and at both temperatures.] 


Material 

Heating (80° F (26.5° C) to 
2000° F (1093.5° C)) 

Cooling (2000° F (1093.5° 
80° F (26.5° C)) 

C) to 


Modulus of 
elasticity, 

E 80°F 

Coefficient of linear 
thermal expansion, 

a 80°-2000°F 

Modulus of 
elasticity, 

E 2000°F 

Coefficient of linear 
thermal expansion, 

a 80°-2000°F 


psi 

GN/m 2 

in./in./°F 

m/m/°C 

psi 

GN/m 2 

in./in./°F 

m/m/°C 

Tungsten a 

50xl0 6 

345 

3xl0~ 6 

5. 4xl0' 6 

37xl0 6 

255 

3xl0~ 6 

5.4xl0~ 6 

80Ni+20Cr b 

30 

307 

C 9. 5 

c 17 

CL 

CO 

d 56 

c 9. 5 

c 17 


a Based on data from ref. 10. 
b Based on data from ref. 11. 
c Based on data from ref. 12. 
d Based on data from ref. 13. 
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TABLE H. - MEAN COEFFICIENT OF LINEAR THERMAL 


EXPANSION FOR SOME REFRACTORY AND 
SUPERALLOY MATERIALS 


[All data obtained from ref. 13, unless otherwise noted.] 


Material 

Coefficient of linear thermal expansion, 
a 80°-2000°F 

in. /in. /°F 

m/m/°C 

Tungsten 

3xl0' 6 

5.4xlO" 6 

Tantalum 

3 

5.4 

Molybdenum 

3.4 

6.1 

80Ni+20Cr a 

9.5 

17 

Rene 41 

10 

18 

Inconel 600 

10 

18 

Haste lloy X 

9.5 

17 

Nimonic 80^ 

10 

18 

GMR 235 

9.8 

17.2 

Nichrotung 

9.5 

17 

L-605 

10 

18 

S-816 

9.2 

16.6 

Wi-52 

9.8 

17.2 

Mar M302 

9.0 

16.2 


^Data from ref. 12. 

^For room temperature to 1600° F (877° C). 



TABLE III. - YIELD AND TENSILE STRENGTHS FOR SOME BULK 


REFRACTORY AND SUPERALLOY MATERIALS 


[All data were based on ref. 13 unless otherwise noted. ] 


Materials 

Temperature, °F (°C) 


80 (26.5) 

2000 (1093.5) 


Strength 


Yield 

Ultimate tensile 

Yield 

Ultimate tensile 


ksi 

MN/m 2 

ksi 

MN/m 2 

ksi 

MN/m 2 

ksi 

MN/m 2 

Tungsten 

— 

— 

160 

1110 

— 

— 

60 

414 

Tantalum 

70 

484 

155 

1070 

10 

69 

20 

138 

Molybdenum 

100 

689 

170 

1170 

40 

267 

80 

551 

80Ni+20Cr 

a 70 

484 

a 117 

a 807 

5 

35 

5 

55 

Rene 41 

40 

276 

180 

1240 

30 

207 

70 

484 

Inconel 600 

35 

242 

124 

900 

5 

35 

8 

55 

Hastelioy X 

60 

415 

110 

780 

— - 

— - 

30 

207 

Nimonic 80 

90 

624 

155 

1070 

5 

35 

5 

35 

GMR 235 

95 

656 

155 

1070 

35 

242 

80 

551 

Nichrotung 

120 

830 

130 

898 

22 

152 

25 

173 

L-605 

80 

551 

160 

1100 

35 

245 

40 

276 

S-816 

80 

551 

140 

918 

— 

— 

10 

69 

WI-52 

65 

449 

90 

624 

18 

125 

20 

138 

Mar M302 

90 

624 

140 

968 

20 

138 

20 



138 


a Based on data from ref. 11. 


TABLE IV. - ELASTIC MODULI FOR SOME REFRACTORY AND 
SUPERALLOY MATERIALS 


[All data based on ref. 13. ] 


Material 

Temperature, °F (°C) 

Material 

Temperature, °F (°C) 


80 (26.5) 

2000 (1093.5) 


80 (26.5) 

2000 (1093.5) 



Elastic moduli 

, E 



Elastic moduli 

E 


ksi 

MN/m 2 

ksi 

MN/m 2 


ksi 

MN/m 2 

ksi 

MN/m 2 

Tungsten 

50 

345 

37 

255 

Hastelioy 

30 

207 

16 

110 

Tantalum 

28 

193 

25 

173 

Nimonic 80 

27 

186 

7 

48 

Molybdenum 

45 

311 

28 

193 

GMR 235 

30 

207 

18 

124 

80Ni+20Cr 

a 30 

207 

8 

55 

Nichrotung 

b 34 

235 

b 22 

152 

Rene 41 

32 

221 

16 

110 

L-605 

34 

235 

20 

138 

Inconel 600 

30 

207 

8 

55 

S-816 

35.5 

245 

24 

166 


a Based on data from ref. 11. 
^Dynamic modulus. 



Tensile and Compressive Stresses on Assumption of Elastic Strain 


Cooling . - Thermally induced equivalent elastic stresses* are presented in figure 

5(a) for a stress-free specimen having been cooled from 2000° F (1093.5° C) to 80° F 

(26.5° C) and in figure 5(b) for a stress-free specimen at 80° F (26.5° C) having been 

heated to 2000° F (1093.5° C). In figure 5(a), it is seen that matrix tensile stresses as 

high as approximately 550 ksi (3790 MN/nr ) can be developed in a high volume-fraction 

reinforcement composite and that reinforcement compressive stresses as high as about 
o 

900 ksi (6200 MN/m ) can be generated in a high volume-fraction matrix composite. 
Tensile and compressive stresses of such high magnitudes could not be achieved because 



.8 .6 .4 .2 0 1.0 

Reinforcement volume fraction, V r 

fa) Cooled from 2000° to 800° F (1093° to 26. 5° C). 


(b) Heated from 80'" to 2000' F. 


Figure 5. - Equivalent elastic stress calculated for model tungsten - 80Ni+20Cr composite. 


*The values of stress and strain must obviously become zero when the volume frac- 
tion of either constituent is zero; that is, the equations are valid only over the interval 
0 < V m r < 1 and must be zero at these limits. 
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they can far exceed both the yield and ultimate tensile strengths of the constituents. 
Yielding or fracturing of either constituent would be expected when its stress became 
sufficiently high. 

Heating. - Upon heating from 80° F(26.5° C)to 2000° F (1093. 5° C). eauivalent elas- 
tic matrix compressive stresses as high as about 150 ksi (1030 MN/m ) could be gener- 
ated for a high volume-fraction reinforcement composite and equivalent elastic reinforce- 
ment tensile stresses as high as about 700 ksi (4820 MN/m 2 ) will be generated in a low 
volume-fraction reinforcement composite. Again, in either instance, yielding or frac- 
turing could occur in the constituent whose yield or ultimate strength has been exceeded. 


Tensile and Compressive Stresses and Strains on Assumption of 
Elastic-Plastic Deformation 

The method and assumptions used in obtaining these values are discussed in appen- 
dix A. 

Cooling . - Thermally induced equivalent elastic-plastic stresses are presented in 
figure 6(a) for specimens having been cooled from 2000° F (1093.5° C) and in figure 6(b) 



1.0 .8 .6 .4 .2 0 1.0 .8 .6 .4 .2 0 

Reinforcement volume fraction, V m 

(at Cooled from 2000° to 80° F (1093. 5° to 26. 5° C). (b) Heated from 80° to 2000° F <26. 5° to 1093. 5° C). 

Figure 6. - Equivalent elastic-plastic stress calculated for model tungsten - 80Ni+20Cr composite. 
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Strain, « 



(c) Tungsten at 2000° F (d) 80Ni+20Cr at 2000° F 

(1093. 5° C) (ret. 10). (1093° C) (ref. 13). 


Figure 7. - Uniaxial stress as function of strain. 



II 


for specimens having been heated to 2000° F (1093. 5° C). The values shown in these 

curves were obtained using a graphical procedure described in appendix A and using the 

stress-strain diagrams of figures 7(a) and (b), which are based on references 10 and 11. 

(A sample calculation is shown in table V. ) In either instance, the stresses, for given 

relative amounts of reinforcement, are less than the calculated elastic stresses. Upon 

o 

cooling, the maximum matrix tensile stress can be as high as about 76 ksi (524 MN/m ), 
and the maximu m reinforcement compressive stress can be as high or about 500 ksi 
(3440 MN/m ) for high volume-fraction reinforcement and high volume-fraction matrix 
composites, respectively. 


TABLE V. - SAMPLE CALCULATION OF ELASTIC-PLASTIC STRESS AND STRAIN 


FOR TUNGSTEN - 80Ni+20Cr COMPOSITE 


[Stresses and strains due to cooling from 2000° F (1093. 5° C) to 80° F (26. 5° C); matrix 
volume fraction, 0.1; reinforcement volume fraction, 0.9; equivalent reinforcement 


plastic strain, e er p is nil for this example. | 


Equivalent matrix stress 

Equivalent 

Equivalent reinforcement stress 

Total 



0. 1 


matrix 



0.9 


equivalent 

a 

em 

a 

em 

plastic 

a 

er 

a er 

plastic 

ksi 

MN/m 2 

ksi 

MN/m 2 

strain, 

e emp 

ksi 

MN/m 2 

ksi 

MN/m 2 

Strain, 

e ept 

110 

758 

11 

78 

0.22 

-12 

-83 

-11 

-7.8 

0.22 

100 

689 

10 

60 

. 14 

-10 

-69 

-9.9 

-6.9 

. 14 

90 

624 

9 

62 

.08 

-9.9 

-69 

-9 

-6.2 

.08 

80 

551 

8 

55 

.04 

-8. 8 

-61 

-8 

-5.5 

.04 

a 78. 5 

541 

7.85 

54 

.0250 

-8.7 

-60 

-7. 85 

-5.4 

.025 

70 

482 

7 

48 

. 002 

-7. 7 

-53 

-7 

-4. 8 

.002 


Value where P„ 


P = 0 
er 


and 


e - e ~ 0. 025 = e , . 
emp erp cpt 


Heating. - Upon heating stress-free specimens from room temperature, the maxi- 
mum matrix equivalent compressive stress is about 7 ksi (48 MN/m ) for high volume- 
fraction reinforcement composites, and the maximum reinforcement equivalent tensile 
stress is about 40 ksi (276 MN/m ) for high volume-fraction matrix composites. The 
stresses generated by cooling or heating would not be expected to cause failure of either 
constituent. However, when heating and cooling cycles are applied, the effect could be 
tantamount to thermal fatigue, as will be pointed out subsequently. 
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Strains and Thermal Fatigue 


Figure 8(a) indicates the approximate magnitudes of equivalent plastic strain which 
occurs when the model laminate composite is cooled from 2000° to 80° F (1093.5° to 
26.5° C), and figure 8(b) indicates the approximate magnitudes of plastic strain when the 
model laminate composite is heated from 80° to 2000° F (26. 5° to 1093. 5° C). Values 
of plastic strain are plotted as a function of the volume fraction of matrix. When values 
become small, plastic strain is plotted as dashed lines to indicate that elastic strain val- 
ues may be significant. Maximum plastic strains as great as about 0.025 have been ob- 
tained for either constituent when the composite is cooled or heated from an initial 
stress-free state. A complete heat-cool cycle can therefore produce a strain range as 
large as about 0.025 in either constituent. At volume fractions of matrix V m up to 
about 0. 80, the strain is essentially occurring in the matrix constituent. With higher 
values of V m O>0. 85), the relative amount of reinforcement decreases to the point 
where the unit stress in the reinforcement is sufficiently high to cause essentially all 
the strain to take place in the reinforcement. Substantial strain will occur in either the 
matrix or reinforcement; if it does not occur in the one, it can occur in the other. 

Data published in references 14 to 16 indicate that steels and alloys can fail in about 
1000 cycles when exposed to a thermally induced total strain range of 0.02. Values of 
strain shown in figure 8 for the model laminate composite produce plastic strain ranges 
that are of the magnitude just noted; hence material failure may occur in about 1000 



1.0 .8 .6 .4 .2 0 1.0 .8 .6 .4 .2 0 

Reinforcement volume fraction, V r 

(a) Cooled from 2000° to 80° F (1093. 5° to 26. 5° 0. (b) Heated from 80° to 2000° F (26. 5 to 1093. 5° C). 


Figure 8. - Equivalent plastic strains calculated for model tungsten - 80Ni+20Cr composite. 
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thermal cycles. The actual number of cycles to failure depends, of course, on test con- 
ditions, including material, the relative amounts of elastic and plastic strain, maximum 
cycle temperature, dwell time at temperature, and, certainly, total strain. 

Of cogent interest are the results obtained in an experiment conducted on S-816 alloy 
gas turbine blades run in a J-47 engine (ref. 17). These blades developed leading-edge 
cracks after 85 normal start-stop cycles (i. e. , heat-cool cycles), involving a thermal 
gradient of 840° F (450° C) at the leading-edge of the blades. Because of this tempera- 
ture gradient, the hotter leading edge was constrained from expanding and was compres- 
sively strained during heating. The reverse took place upon cooling. The number of 
cycles to initiate cracking was lessened by introducing periods of normal operation. In 
a subsequent paper (ref. 18) the elastic strain imposed upon the leading edge of these 
blades, was calculated to be somewhat less than 0.004. Thus a total strain range of 
about 0. 004 may have been imposed on these blades during a start-stop cycle. 

Laboratory test data indicate relatively rapid thermal fatigue damage (i. e. , failure 
in 1000 cycles) for numerous alloys subjected to strain ranges on the order of 0.02. 
However, in an actual application, a calculated elastic strain range of about 0.004 was 
associated with thermal fatigue cracking in 85 heat-cool cycles; it may be speculated 
that creep (i. e. , dwell time under stress) may have been a contributing factor in causing 
the aforementioned cracking. The preceding data and conjecture very strongly suggest 
that, when the total strain range approaches a magnitude of 0.02 per cycle, thermal fa- 
tigue failure may ensue in relatively few cycles (i. e. , 1000). The question may also be 
asked whether or not any losses in mechanical properties may occur prior to any visible 
material failure; this may be particularly important in regard to the reinforcement ma- 
terial. Thus the use of sheet- or foil- reinforced composites (or fiber- reinforced com- 
posites with multidirectional plies) might pose thermal fatigue problems since, in actual 
usage where thermal cycling might be involved, the strain ranges encountered may be a 
magnitude sufficiently high to cause material or mechanical property deterioration. 

(The thermal expansion of a fiber-reinforced ply or tape can depend on direction, 
being least in the direction of the fiber if the fibers have the lower coefficient of expan- 
sion and being greatest in the direction orthoganal to that of the fibers. Therefore, two 
fiber- reinforced plies, joined to each other so that the fibers in one ply were at some 
angle to the fibers in the other ply, would be analogous to alternate sheet or foil of dif- 
fering thermal expansivities. The type of analysis derived in this study for sheet and 
foil composites would then be applicable to this particular type of composite insofar as 
average stresses and strains are concerned. ) 
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Elastic Shear Stress 


The calculated peak elastic shear stresses resulting from the cooling or heating of 
the stress-free model laminate composite are shown in figure 9. These values are shown 
as a function of the volume fraction of constituent. These stresses occur very near 
free edges and diminish as one moves away from the edge. The calculated stresses are 
conservative when V <0.5 (see appendix F). On cooling, maximum shear stresses of 

about 140 ksi (965 MN/m^) occur at V ~ 0. 7; over the range of 0. 2 < V <0.9, the 

shear stresses are above 100 ksi (689 MN/m ). On heating, the maximum peak shear 

o 

stresses approach 40 ksi (276 MN/m ) at V_ « 0. 8 and are in excess of 20 ksi (138 
n m. 

MN/in ) over the range of 0. 2 < V m< 0.9. Assuming the shear strength to be one-half 
the tensile strength, the shear strengths may be estimated to be 60 and 2.5 ksi (415 and 
17.3 MN/m^) for the 80Ni+20Cr alloy at 80° and 2000° F (26.5° and 1093.5° C), respec- 



1.0 .8 .6 .4 .2 0 

Matrix volume fraction, V m 

Figure 9. - Peak elastic shear stresses calculated for matrix of tungsten - 
80Ni+20Cr laminate composite after heating or cooling between 80° and 
2000° F (26. 5° and 1093. 5° C). 
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tively. These estimates of shear strength are based on the tensile values shown in fig- 
ure 7. Thus, it is possible that shear failure could occur in the matrix (or at the inter- 
face) particularly over the range of volume fractions of matrix that might be used in 
practice. 

An observation of potential importance may be made. From an inspection of table 

III, it may be seen that the ultimate tensile strengths of a number of candidate matrix 

o 

materials is a maximum of about 160 ksi (1104 MN/m ) at room temperature and about 
80 ksi (551 MN/m^) at 2000° F (1093. 5° C). If shear strengths are estimated to be one- 
half the tensile strengths, shear strengths for matrix materials may approximate 80 ksi 
(551 MN/m^) at room temperature and 40 ksi (276 MN/rn^) at 2000° F (1093. 5° C). Fur- 
thermore, the elastic moduli of many potential matrix materials (see table IV) are simi- 
lar to or greater than that of 80Ni+20Cr so that the shear moduli would be expected to be 
similar or greater. Thus, there is a likelihood that the peak shear stresses in many 
possible matrix materials will exceed their shear strengths. 

If plastic deformation occurs rather than fracture, the stresses in the matrix will 
be less than calculated, but the effect of yielding (particularly cyclic yielding) on bond 
integrity is not known and repeated yielding of the matrix due to shear stress may cause 
bond separation or degradation. 

A very important point can be made; namely, if discontinuities are introduced into 
the composite (i. e. , if the composite is altered by machining or any other process), free 
edges may be introduced. Because the highest thermally induced shear stresses occur 
at or near free edges, these modifications can introduce more potential shear failure 
sites into the material. Thus, from this particular standpoint, changes in cross section, 
cutouts, contours, reinforcement or matrix breaks, porosity in the matrix, etc., are all 
potentially deleterious . 


EXPERIMENTAL STUDY OF LAMINATES 

A study of limited scope was conducted to obtain indications of whether readily ap- 
parent damage could be induced in reinforced sheet or foil laminate composites. The 
calculated stresses in the tungsten - 80Ni+20Cr composite model were determined to be 
sufficiently high to result in tensile, compressive, or shear failure during the first cool- 
ing or heating cycle if plastic stress relief did not occur. And should stress relief occur, 
thermal fatigue might later ensue after a number of heat-cool cycles. Simple heat-cool 
tests were used to determine whether a relatively few heat-cool cycles could produce 
readily apparent laminate composite deterioration. There was no intent to necessarily 
induce thermal fatigue in the event this type of damage should require large numbers of 
cycles. Heat-cool cycles with very slow heating and cooling rates were used, as well as 
heat-cool cycles with fairly rapid heating and cooling rates. Thermal gradients would be 
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minimal for slow heating and cooling rates, if they occurred at all. But in actual use con- 
ditions, such as might be encountered in a turbine blade application, rapid heating and 
cooling could produce temperature gradients and, therefore, produce stresses which are 
additive to those due to the differences in the thermal expansions of the constituents. 


Materials 


Commercial tungsten and 80Ni+20Cr sheets and foils were used as the laminae. 

Table VI shows the thicknesses of the laminae used. The specimens were formed by 
alternating layers (lamina) of tungsten and 80Ni+20Cr. Six specimens approximately 1- 
by 3-inch (2. 54- by 7. 12-cm) thick were tested. Three were as-fabricated (rectangular 
specimens, fig. 10(a)), and three had the configuration shown in figure 10(b) (contoured 
specimens). The rectangular specimens were used in the slow heat-cool tests; the con- 
toured specimens in the more rapid heat-cool tests. Both types of specimen had 0.001-, 
0.005-, and 0. 020-inch (0. 0025- , 0.0125-, and 0. 050-cm ) thick tungsten laminae with cor- 
responding 80Ni+20Cr laminae thicknesses (see table VI). The tungsten laminae were 
used in the as-received thickness; the 80Ni+20Cr were obtained in 0.001- and 0.005-inch 
(0.0025- and 0.0125-cm) thicknesses, and laminae over 0.005 inch (0.0125 cm) thick 
were made by using multiple layers of 0.005-inch (0.0125-cm) thick foils. 

The volume fraction of tungsten was about 0.50 in all instances. The different foil 
thicknesses were used because they were conveniently available and because any possible 

TABLE VI. - SPECIMEN DETAILS 


[All thickness values are approximate.] 


Specimen shape 

Tungsten laminae 

80Ni+20Cr laminae 

Total 

number 

of 

Specimen thickness 


Thickness 

Number 

Thickness 

Number 

in. 

cm 

Rectangular a 

in. 

0.001 

cm 

0.0025 

20 

in. 

0.001 

cm 

0.0025 

21 

laminae 

41 

0.041 

0.1 


.005 

.0125 

10 

.005 

.0125 

11 

21 

. 105 

.263 


.020 

.050 

4 

b . 020 

b .050 

5 

9 

. 180 

.457 

Contour 0 

0.001 

0.0025 

20 

0.001 

0.0025 

21 

41 

0.041 

0.1 


.005 

.0125 

10 

.005 

.0125 

11 

21 

. 105 

.263 


.020 

.050 

2 

(d) 

(d) 

3 

5 

.080 

.200 


a Used for slow heat- cool cycles. See fig. 13(a). 

b Four 0. 005-in. (0. 0125-cm) thick 80Ni+20Cr foils were used to form the 0. 020-in. (0. 0125-cm) 


thick lamina. 

c Used for rapid heat-cool cycles. See fig. 13(b). 

^The outer 80Ni+20Cr laminae were 0.010 in. (0.025 cm) thick (two 0.005-in. (0.0125 cm) foils); 
the inner lamina was 0.020 in. (0.050 cm) thick (four 0.005-in. (0.0125-cm) foils). 
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(a) Rectangular specimen used for slow heat-cool tests. 



Figure 10. - Configurations of heat-cool test specimens. 
(Dimensions are in inches (cm).) 


effect of lamina thickness on composite susceptibility to damage was of interest. All 
specimens were consolidated (bonded) by vacuum hot pressing for 4 hours at 1800° F 
(980° C) and 4 ksi (28 MN/m^) or better. Specific details of the specimens are given in 
table VI. 


Procedure 

The specimens having the rectangular profile were exposed simultaneously (i. e. , 

three at a time); they were slowly heated in a vacuum furnace for 2 hours to 2000° F 

(1093.5° C) and allowed to remain at this temperature for 2 hours. The specimens 

were then slowly cooled to room temperature for 2 hours. The vacuum was 10 - ^ torr 
o 

(13 MN/m ) or better. They were inspected after each heat-cool cycle, using magnifi- 
cations up to X150. 

The contoured specimens (fig. 10(b)) were each placed into a separate heat treating 
envelope. Each specimen was placed in a furnace which had been preheated to 1600° F 
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(877° C) and allowed to remain in the furnace for 15 minutes. Each specimen was given 
its first cooling cycle by being immersed in a container of water. But the envelopes be- 
came filled with water, and subsequent cooling was accomplished by placing the specimen 
on a metal plate which was at ambient temperature. The specimen having the 0.020- 
inch tungsten laminae failed after the first cooling cycle and was not tested further. The 
160 0° F (877° C) furnace temperature was used because of the furnace temperature limi- 
tations and because it was desired to minimize specimen oxidation which might occur. 

Varying numbers of cycles were used in both the slow and fast cycle tests in order to 
obtain a balance between observable effects and ready availability of evaluation apparatus 


Results and Discussion 

Slow heat-cool cy cles. - The results of the slow heat-cool cycles are summarized 
in table VII(a). The composite with the 0.020- (0.050-cm) thick laminae, exhibited de- 
lamination after the first cycle and progressive delamination up to the third cycle. The 
damage appeared to be essentially stabilized at this time. The specimen with the 0. 005- 
inch (0.0125-cm) laminae exhibited microdamage after 11 heat-cool cycles. The speci- 
men with the 0.001-inch (0.0025-cm) laminae did not exhibit any observable damage 
through the 11th heat-cool cycle. The appearance of the specimens after 11 cycles is 
shown in figure 11. 

Even with very slow heating and cooling, visible damage did occur. The calculated 

o 

elastic-plastic stresses for the matrix were 7 and 76 ksi (48. 2 and 523 MN/m ) on heat- 
ing and cooling, respectively, assuming that relaxation did not occur. The calculated 
stresses in the reinforcement were also 7 and 76 ksi (48. 2 and 523 MN/m ) on heating 
and cooling, respectively. The calculated plastic strain range in the matrix was about 

0. 025 per heat-cool cycle. The calculated elastic shear stress in the matrix was 30 and 

2 

140 ksi (208 and 965 MN/m ) during heating and cooling, respectively. With a strain 
range of about 0. 025, thermal fatigue damage might be expected in less than 1000 cycles, 
and, possibly, thermal fatigue damage actually occurred in the specimen with the 0. 005- 
inch (0. 0125-cm) laminae. 

Rapid heat-cool cycles . - The results of this experimentation are summarized in 
table VII (b). The specimen containing the 0.020-inch (0.050-cm) tungsten laminae was 
completely unbonded after the first cycle. The specimen containing the 0. 005-inch 
(0.0125-cm) tungsten laminae, which showed no damage after the first heat-cool cycle, 
was very badly damaged after two more cycles. The laminate specimen containing the 
0.001-inch (0.0025-cm) tungsten laminae exhibited incipient laminae separation after 
six heat-cool cycles. The appearance of each of these specimens at the end of its re- 
spective final test cycle is shown in figure 12. The calculated laminae stresses and 
strains correspond to those previously indicated for the slow-cool specimens. Because 
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TABLE VII. - RESULTS OF FURNACE THERMAL CYCLING 


(a) Slow cool 


Laminae thickness 

Number 

of 

cycles 

Results 

(a) 

Tungsten 

80Ni+20Cr 

in. 

0.001 

cm 

0.0025 

in. 

0.001 

cm 

0.0025 

11 

No macroscopically visible damage 

0.005 

0.012 

0.005 

0.012 

11 

Cracking of lamallae - possible delamination 

0.020 

0.050 

0.020 

0.050 

1 

Delamination 





2 

Further delamination 





3 

Further delamination 





11 

Slightly further delamination 


(b) Rapid cool* 3 


Laminae thickness 

Number 

Results 

Tungsten 

80Ni+20Cr 

of 

cycles 


in. 

cm 

in. 

cm 


(c) 

0.001 

0.0025 

0.001 

0.0025 

6 

Incipient foil separation 

0.005 

0.0125 

0.005 

0.0125 

3 

Foil separation, shear buckling 

0.020 

0.050 

(d) 

(d) 

1 

Sheets unbonded 


a See fig. 11. 

^First quench was into water container. Subsequent quenches by air cool on metallic 
plate at ambient temperature. Specimens were in a 1600° F (877° C) furnace for 
15 min and then removed for cooling. 
c See fig. 12. 

d See footnote (d), table VI. 
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(a) Tungsten lamina thickness, 0.020 inch (0.050 cm). Appearance essentially 
unchanged from that observed after three cycles. Note gross delamination of 
outer elements and separation of inner elements. 



(b) Tungsten lamina thickness, 0.005 inch (0.0125 cm). Note cracking of both (c) Tungsten lamina thickness, 0.001 inch (0.0025 cm). No discernable 

tungsten and 80Ni+20Cr laminae. There is indication of bond separation. damage. 

Figure 11. - Tungsten - 80Ni+20Cr composite after 11 slow heat-cool cycles to and from 80° and 2000° F (26.5° and 1093.5° 0. Tested in as-hot-pressed condition. 
Photographs reduced by 20% in reproduction. 



(a) Tungsten lamina thickness, 0.020 inch (0.050 cm); after one cycle. 


(b) Tungsten lamina thickness, 0.005 inch (0.0125 cm); after three cycles. Note shear failure at 
interface, breaking, and damage to tungsten lamina. 



Figure 12. - Tungsten - 80Ni+20Cr composite before (topiand after (bottom) rapid heat-cool cycles 
to and from room temperature and 1600° F (877° 0. Specimen tested in as-hot-pressed condi- 
tion. X2. 









KP"' • . 

j, \ 

''—Incipient damage 

0-70-1123 

(C) Tungsten lamina thickness, 0.001 inch (0.0025 cm); after six cycles. 

Figure 12. - Concluded. 

of the relative rapidity of temperature change, there is a possibility of introducing ther- 
mal gradients and concommitant additive stresses, so that the end result in these par- 
ticular tests are not necessarily attributable to thermal expansion differences alone. 
Actually, however, rapid temperature changes could well be encountered in practical 
applications, so that these particular results obtained in the relatively rapid heat-cool 
study could indicate what might actually be encountered. 

Regardless of whether the heating and cooling was slow or rapid, the greatest ob- 
served structural damage occurred to the specimens with the thickest laminae, and the 
least to the specimens with the thinnest laminae. In both instances the 0. 020-inch 
(0. 050-crn) laminae specimen sustained delamination, which indicates a predominance of 
shear effects. Inspection of the equation for shear stress (eq. (3) or (B21)) indicates that 
the shear stress decreases more rapidly away from the tip, as the matrix laminae thick- 
ness decreases. Thus, the region of high shear stress would be greatest for the 0.020- 
inch (0.050-cm) thick laminae composites and least for the 0.001-inch (0.0025-em) thick 
laminae composites; and it might explain the apparent decreased effect of shear with de- 
creased laminae thickness. Another possible reason for lesser damage in the thinner 
laminae composites is that a given amount of interdiffusion would involve a relatively 
larger amount of each lamina. This could tend to lessen the differences between the coeffi- 
cients of thermal expansion of adjacent lamina and thereby lessen the stresses actually 
developed. This effect would be greatest forthe 0. 001-inch (0.0025-cm) laminae compos- 
ites. Obviously, the more rapid temperature change accelerated the structural damages. 


SIGNIFICANCE OF MATHEMATICAL ANALYSIS AND EXPERIMENTAL RESULTS 


The mathematical stress analysis and the resulting calculated stresses for the 
tungsten - 80Ni+20Cr laminate composite indicated that very high tensile, compressive, 
or shear stresses (in excess of the ultimate tensile and the estimated shear strengths of 
the materials) could be generated because of the difference in the coefficients of thermal 
expansion of the constituents and because of a temperature transition between 2000° and 
80° F (1093.5° and 26.5° C). The materials combination and the temperature ranges 
were representative of a genre of composite materials and use- conditions that are being 
considered for practical application. The high calculated elastic stresses could be me- 
liorated by plastic deformation but the concurrent strains might cause thermal fatigue 
damage if cyclic heating and cooling were encountered. From an analytical viewpoint, 
the model system studied seemed prone either to immediate structural damage due to 
high stresses or to limitations on its endurance due to strain associated with cyclic ther- 
mal stresses. 

An experiment involving relatively few heat- cool cycles applied to tungsten - 
80 Ni+20Cr laminate composites resulted in damage (whether the heating and cooling 
rates be slow or rapid) for laminae thicknesses that might be of practical interest. The 
extent of this damage was greater with increased laminae thicknesses and with more 
rapid heating and cooling rates. 

Thus, sheet or foil laminate composites may be prone to damage as the result of 
high thermally induced stresses arising from differences in coefficients of thermal ex- 
pansion of the constituents and exposure to one or more temperature transitions. 


FIBER COMPOSITES 

STRESSES IN FULL-LENGTH FIBER-REINFORCED COMPOSITES 

Fiber-reinforced composites have been arbitrarily catagorized as being high volume 
fraction if they have between 65 and 90 percent fiber and low volume fraction if they have 
less than 65 percent fiber content. The 90 percent figure corresponds to the theoretical 
packing limit of circular fibers. The 65-percent figure is about the minimum amount of 
reinforcement assumed to be usually associated with high volume-fraction composites. 

The low volume-fraction fiber composites were assumed to be amenable to analysis 
using thick shell theory, and the high volume-fraction fiber composites were assumed to 
be in a state of hydrostatic stress. A rigorous mathematical analysis incorporating the 
effect of volume fraction of fiber and related constraints was beyond the scope of this 
study. 
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High Volume-Fraction Composites (0.65<V r <0.90) 


Assumptions . - 

(1) The fibers are well bonded to the matrix. 

(2) Each constituent is homogeneous and isotropic, and has the same Poisson's ratio. 

(3) A hydrostatic stress state exists within each constituent of the composite; that 
is, = Or = cr x . The hydrostatic stress state is the result of a sufficiently high volume 
fraction of fiber so that there is stress interaction between matrix and fibers. 

(4) Tensile and compressive stresses are uniform within the breadth of the respec- 
tive constituents. 

(5) The fibers are in a hexagonal array. 

(6) There is no interphase region. 

(7) End effects are neglected. 

Mechanism of stress generation . - A model is shown in figure 13. Figure 13(a) il- 
lustrates a complete composite containing many full-length fibers. Figure 13(b) repre- 
sents an element of this composite comprising a fiber and an "associated” shell of ma- 
trix. The polar coordinate system used is shown in figure 14. The fiber and matrix 
have different coefficients of thermal expansion, and the composite is assumed to have 
been consolidated at elevated temperature and cooled to room temperature. Because of 



(a) Complete composite segment 


°tension' m-* 



(c) Fiber-matrix element from low 
volume fraction fiber composite. 


tension, m-\ 



(d) Cylinder within a showi ng fiber, matrix, and polar coordinate 

cylinder model.. system. 


Figure 13. - Fiber-reinforced composites. 
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the greater thermal contraction of the matrix, it will contract about the fiber. Stresses 
will be introduced into the fiber and the matrix. These stresses will be in the longitudi- 
nal (as shown), radial, and tangential directions. If a consolidated composite is heated, 
stresses of an opposite sign to those produced on cooling will result. Stresses and 
strains will be produced on each subsequent heating or cooling. The derivation of equa- 
tions for longitudinal, radial, and tangential stress is given in appendix C. The resul- 
tant equations are as follows: 


°xm °Rm a 6m 


VE E (a - a )AT 
r r m x r nr 

(V E + V E )(1 - 2v) 
\ » r -^ r T * m m /v ’ 


u = ov» = = 

xr Rr 6x 


V E E (a -a )AT 
m m r v m r' 

(V E + V E )(1 - 2v) 
' r r m m /v ' 


(5) 


( 6 ) 


There will be no plastic deformation since it has been assumed that a hydrostatic stress 
state exists. 


Low Volume-Fraction Composites 


Assumptions . - 

(1) Matrix reinforcement are perfectly bonded. 

(2) Matrix and reinforcement are, respectively, homogeneous and isotropic and have 
the same Poisson's ratio. 

(3) There is no stress interaction between matrix and neighboring fibers; thus, the 
thick-wall cylinder theory applies. 

(4) The fibers are in a hexagonal array. 

(5) There is no interphase region. 

(6) End effects are neglected. 

Mechanism of stress generation . - The model is represented in figure 13, an ele- 
ment of which is illustrated in figure 13(c). The coordinate system is shown in figure 14. 
The composite, consisting of matrix and fibers with different coefficients of thermal ex- 
pansion is assumed to have been consolidated at elevated temperature and cooled to room 
temperature. Both constituents are initially the same length and after cooling are again 
of equal lengths. This will generate longitudinal (as shown), radial, and tangential 
stresses in fiber and matrix as before, but, because of assumption (3), plastic flow can 
take place. If the composite were stress free at room temperature and heated, stresses 
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would again be generated, but of opposite signs. The derivations are presented in ap- 
pendix D and the resulting equations are 
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and b and c are fiber and associated matrix radii, respectively, as shown in figure 
13(d). Equations (7) and (8) give stresses at the fiber-matrix interface; the matrix 
stress are theoretically the largest at this point. 

Equations (7) to (12) yield elastic stress values but may be used to obtain estimates 
of elastic-plastic stresses and strains as discussed in appendix C. 
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ELASTIC SHEAR STRESSES IN FULL-LENGTH FIBER-REINFORCED COMPOSITES 


Assumptions 

An equation for shear stress may be derived in a manner analogous to that used for 
the laminate composite. All previous assumptions made for the laminate composite still 
apply with the exceptions or additions that 

(1) A triaxial rather than a biaxial stress-state exists. 

(2) The fibers are in a hexagonal array. 

(3) The interfiber distance is equal to or less than the fiber diameter; this corre- 
sponds to V r > 0.23. 



x 

(a) Typical element (b) Forces on reinforcement element. 


Figure 15. - Element of fiber reinforced composite and forces acting on 
reinforcement (i.e., fiber). 


Mechanism of Stress Generation 


A plane section through the fiber and its associated cylinder of matrix may be repre- 
sented by figure 15. The comments made regarding stress generation (p. 7) apply, and 
the actual derivation is given in appendix E. The resulting equation for elastic stress is 
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CALCULATION OF RESIDUAL THERMAL STRESSES AND STRAINS 
IN MODEL FIBER-REINFORCED COMPOSITE 

Residual stresses have been calculated for a model fiber- reinforced composite con- 
taining tungsten fibers in a 80Ni-20Cr matrix, as was the case for the laminate compos- 
ites. Conditions A and B (p. 9) will be used again as will the respective material proper- 
ties as shown in table I. The purpose of these calculations is to obtain representative 
stress and strain magnitudes for the model composite, and these values will be applica- 
ble to other materials combinations having similar coefficients of thermal expansion, 
Poisson's ratios, yield strengths, and elastic moduli. The results of the calculations 
are shown in figures 16 to 21. 


Tensile and Compressive Stresses in High-Volume-Fraction 
Fiber-Reinforced Composites 

Cooling . - Fiber contents ranging from an arbitrary 0.65 to 0.90 volume fraction 
(the maximum theoretical fiber content for circular fibers is 90.3 percent) have been 
assumed to comprise high volume-fraction fiber- reinforced composites. Stresses for 
composites with V r > 0. 90 have been indicated by dashed lines since V r cannot exceed 
about 0.90 for circular fibers. Figure 16(a) indicates that on cooling, the tensile stress 
(i. e. , Ug, O]^, or ct z ) in the matrix can range from about 750 to about 600 ksi (5170 to 
o 

4150 MN/m ). Tensile stresses of this magnitude should cause fracture of the matrix 

(unless yielding occurs). The calculated compressive stress in the fibers may be as 

2 

high as 350 ksi (2430 MN/m ) (for the 0.65 volume fraction fiber composites), and com- 
pressive stresses of this magnitude should not cause any immediate harm to the tungsten 
reinforcement. 

Heating. - Figure 16(b) indicates that matrix compressive stresses (i. e. , o a , a D , 

2 o K 

and a ) as high as 200 ksi (1378 MN/rn ) may be generated on heating. These may frac- 
ture the matrix (unless relieved by plastic yielding). Tensile stresses as great as 
100 ksi (689 MN/m ) may be generated in the reinforcement. Stresses of this magnitude 
would probably fracture the tungsten fibers (unless plastic deformation occurred). 
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1.0 .8 .6 1.0 
Reinforcement volume fraction, V r 


.6 


(a) After cooling from 
2000° to 80° F 
(1093. 5° to 26.5° 0. 


(b) After heating from 
80° to 2000° F (26. 5° 
to 1093. 5° 0. 


Figure 16. - Axial elastic stresses (o x = og » or) calculated for model high-volume- 
fraction (0. 65 < V r < 0.90) tungsten - 80Ni+20Cr-matrix fiber composites. 


Tensile and Compressive Stresses in Low-Volume-Fraction 
Fiber-Reinforced Composites 

Cooling. - Figure 17(a) indicates that longitudinal compressive stresses as great as 

2 

about 700 ksi (4716 MN/m ) can be generated in the fibers and tensile stresses as great 

2 

as about 250 ksi (1722 MN/m ) could be generated in the matrix. Transverse tensile 

stresses (fig. 17(b)) ranging from about 250 to about 380 ksi (1722 to 2600 MN/m^) could 

be developed in the matrix. Transverse compressive stresses as high as about 250 ksi 
2 

(1722 MN/m ) could be developed in the matrix and reinforcement with low matrix vol- 

o 

ume fractions; these stresses decrease to about 75 ksi (516 MN/m ) with a reinforce- 
ment volume fraction of about 0. 65. The reinforcement stresses might not be damaging, 
but the matrix stresses could cause fracture on first cooldown if not relieved by plastic 
flow. 

Heating. - Figure 18(a) indicates longitudinal fiber tensile stresses as high as about 
2 

450 ksi (3100 MN/m ) at small reinforcement volume fractions; these stress values de- 

2 

crease to about 40 ksi (276 MN/m ) with an increase in the reinforcement volume frac- 
tion to about 0. 65. The compressive stress in the matrix increases to about 80 ksi 
(55 MN/m ) as the reinforcement volume fraction increases to about 0. 65. The trans- 
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1.0 .8 .6 .4 

(a) Longitudinal stress. 


.2 1.0 
Matrix volume fraction, V m 


.8 .6 .4 .2 

(b) Transverse stress calculation for interface. 


Figure 17. - Elastic stresses calculated for model low volume fraction ( V r < 0. 65) tungsten-fiber-reinforced - 80Ni+20Cr-matrix 
composites. After cooling from 2000° to 80° F (1093. 5° to 26. 5° 0. 


verse compressive stresses (fig. 18(b)) in the matrix increase from about 75 ksi (503 
MN/rn ) at low reinforcement volume fractions to about 100 ksi (689 MN/m ) at matrix 
volume fractions of about 0. 65. The reinforcement and matrix tensile stresses decrease 
from about 80 to about 20 ksi (551 to 138 MN/m ) over this same reinforcement volume- 
fraction interval. Microfracture of either constituent could occur upon first heating at 
some of these stress levels unless relieved by plastic deformation. 


Tensile and Compressive Stresses in Low-Volume -Fraction Fiber-Reinforced 
Composites on Assumption of Elastic-Plastic Deformation 

These stresses and strains were calculated using the method and assumptions de- 
scribed in appendix D. 

Cooling. - Equivalent stresses were calculated and are shown in figure 19. A max- 
imum compressive stress in the fiber of about 300 ksi (2070 MN/m ) may be obtained at 
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Matrix volume fraction, V m 


(a) Longitudinal stress. (b) Transverse stress calculation for interface. 

Figure 18. - Elastic stresses calculated for model low volume fraction (V r < 0.65) tungsten-fiber-reinforced - 80Ni+20Cr~matrix 
composites. After heati ng from 80° to 2000° F (26. 5° to 1093. 5° 0. 



Reinforcement volume fraction, V r 



Matrix volume fraction, V m 

(a) Tungsten fibers. (b) 80Ni + 20Cr matrix. 

Figure 19. - Equivalent elastic-plastic stresses calculated for model low volume fraction (V r < 0. 65) tungsten-fiber-reinforced - 
80Ni + 20Cr-matrix composite. After heating or cooling between 80° and 2000° F (26. 5° and 1093. 5° 0. 
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low reinforcement volume fractions. With reinforcement volume fractions of about 0.65, 

2 

the stress decreases to about 50 ksi (345 MN/m ). The matrix tensile stress is about 
o 

75 ksi (518 MN/m ) for the low reinforcement volume-fraction composites. This stress 

o 

increases to about 80 ksi (551 MN/m ) for reinforcement volume fractions approaching 
0.65. 

Heating. - Inspection of figure 19 indicates that the maximum tensile fiber stress 
2 

will be about 35 ksi (220 MN/m ) for low volume-fraction composites decreasing to about 
2 

20 ksi (138 MN/m ) for a fiber fraction of about 0. 65. The matrix compressive stress 

2 

will be about 6 ksi (42 MN/m ) for composites with fiber volume fractions up to 0.65. 

Strains and Thermal Fatigue in Low-Volume-Fraction 
Fiber-Reinforced Composites 

Elastic-plastic equivalent strains are presented in figure 20. The maximum total 
and plastic strains occurring in the fiber on cooling are about 0.010 and 0.005, respec- 
tively, and these values decrease rapidly as the reinforcement volume fraction is in- 
creased. The maximum total and plastic strains occurring in the fiber on heating are 
0.01 and 0.008, respectively, and also decrease rapidly with increasing reinforcement 
volume fraction. The total and plastic strains in the matrix are quite similar to each 
other. On cooling, the matrix plastic strain is about 0. Oil up to a reinforcement volume 
fraction of about 0.50. Above this value of reinforcement, the matrix strain increases 
rapidly to about 0.033 at a reinforcement volume fraction of 0. 65. On heating, the plas- 
tic strain in the matrix is about 0.014 and remains at this level, even up to reinforce- 
ment volume fractions of 0. 65. The strains in the reinforcement may be damaging when 
the amount of the reinforcement is relatively small (i. e. , V r < 0. 25). Cyclic strain, as 
small as 0. 004, which may occur when V r < 0.25, was associated with gas turbine blade 
cracking (refs. 17 and 18) in only 85 thermal cycles. The strains calculated for the ma- 
trix are of even greater concern because of their greater magnitude and because they 
occur for virtually all amounts of fiber content up to V r = 0.65. Further, an apprecia- 
ble amount of data (refs. 14 to 16), indicates that a cyclic strain range of 0. 02 can cause 
thermal fatigue in 1000 cycles. Since the matrix strain can be on the order of 0.01, ma- 
trix endurance could be cycle limited, even if it exceeds 1000 cycles. In actuality, the 
calculated strains would probably be approached with lower volume fractions of rein- 
forcement, but, as the amount of reinforcement is increased, the amount of strain would 
become progressively lower, until, at high values of reinforcement volume fraction, 
there would be essentially complete restraint and structural damage would occur rather 
than plastic flow. 

Obviously, a refractory- metal-fiber-reinforced alloy matrix composite can be 
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(a-2) Cooling. 


!b-2) Heating. 

(a) Tungsten fibers. 


(b) 80Ni+20Cr matrix. 


Figure 20. - Equivalent total and plastic strains estimated for model low volume fraction (V r < 0. 65) 
tungsten-fiber-reinforced - 80Ni+20Cr-matrix composite. After cooling between 80° and 
2000° F (26. 5° and 1093. 5° C). 


harmed by cyclic thermal stressing. Possible loss of fiber strength or even fracture 
could occur after numbers of thermal fatigue cycles, while thermal fatigue cracks in the 
matrix would permit oxidation of the reinforcement. 


Elastic Shear Stresses 

In order to evaluate equation (15) for shear stress, it is necessary to determine the 
appropriate values of K where 


K 


(a 0 + a R ) 


(17) 


By referring to figures 17 and 18, it is possible to select corresponding values of a 0 , 
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ov,, and a for given values of V and to determine values of K for V < 0.65. For 

V r > 0. 65, K = 2 because of the assumption of hydrostatic stress. Values of peak shear 

stress in the matrix of the model fiber- reinforced composite are shown in figure 21. 

The maximum values occur at a fiber fraction of approximately 0. 80 both on heating and 

o 

cooling. This maximum is about 135 ksi (935 MN/m ) on cooling and about 40 ksi (276 
2 

MN/m ) on heating. It is apparent that the elastic shear stresses in the matrix are 

2 

equal to or in excess of the estimated matrix shear strengths (60 ksi (413 MN/m ) at 
80° F (26. 5° C) and 25 ksi (172 MN/m 2 ) at 2000° F (1093. 5° C)) essentially over the 
practical range of fiber contents. The stresses for composites with a fiber content be- 
low about 23 percent will be conservative for reasons explained in appendix F. 

Incidentally, because A is greater for the fiber composite than for the laminate 
composite, the tip shear stresses in the fiber composite will decay more rapidly away 
from the tip than in the instance of the laminate composite. 

All of the discussion pertaining to the effects of shear stress in laminate composites 
is, of course, germane to fiber-reinforced composites. 



1 I __L L . I I 

1.0 .8 .6 .4 . 2 0 

Matrix volume fraction, V m 

Figure 21. - Peak elastic shear stresses in tungsten-fiber-reinforced - 
80Ni+20Cr-matrix composite. After heating or cooling between 80° 
and 2000° F (26. 5° and 1093. 5° C). 
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EXPERIMENTAL STUDIES OF FIBER-REINFORCED COMPOSITES 


A study of limited scope was run on fiber- reinforced metallic-matrix composites to 
determine whether the fiber- reinforced composites could be damaged by heat-cool cy- 
cles; that is, to determine whether cyclic variations of temperature, which an actual 
fiber- reinforced composite might be exposed to, would produce microstructural evi- 
dence of composite damage. As before, the concern was with possible gross effects 
that might become evident in a limited number of cycles rather than a concern with rela- 
tively large numbers of cycles that might be required to produce thermal fatigue. 


Materials, Apparatus, and Procedure 

Two specimens of tungsten fiber- reinforced copper were evaluated; one specimen 
contained about 0.40 volume fraction tungsten and the other contained about 0. 70 volume 
fraction tungsten. These particular composites were used because they were available 
from a previous study conducted at this laboratory and because they were felt to be met- 
allurgically sound specimens. Further, the coefficient of thermal expansion of copper 
approximates that of nickel base alloys; the copper is presumably quite ductile and is 
theoretically nonreactive with tungsten. The last two characteristics would suggest that 
this system would tend to minimize any possible cyclic effects. The fiber was nominally 
0. 005 inch (0. 0125 cm) in diameter. The composites were made by an infiltration tech- 
nique such as described in reference 19. The specimens were initially about 1-inch 
(2.5-cm) long by 1/3 inch (0. 13 cm) in diameter. The 0. 40- volume-fraction tungsten 
composite had macroporosity. The 0. 70-volume-fraction tungsten composite appeared 
free of macroscopic porosity. A small disk was taken from each of the two composite 
specimens. The larger portion of each composite was polished on either end to aid in 
subsequent examinations. Each disk was metallographically polished on one of its sur- 
faces. 

Each of the larger specimens was placed into a metallic foil heat-treating envelope 
along with a small amount of aluminum oxide (AlgOg) powder to prevent adherence of the 
specimen to the envelope. Each envelope was folded and then inserted into second en- 
velopes. A thermocouple was inserted between the envelopes of one of the specimens to 
monitor temperature. The specimens were inserted simultaneously into an air atmo- 
sphere furnace which had been preheated to 1600° F (877° C). This temperature was 
purposely chosen to keep the test temperature well below the melting point of the copper 
matrix. The total time in the furnace for each heating cycle was 20 minutes, of whF’h 
time 18 minutes were required for the thermocouple to reach 1600° F (877° C). The 
specimens, still in their protective envelopes, were cooled to ambient temperature; 
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they were removed from the furnace and placed on refractory brick. A temperature of 
180° F (76° C) was indicated after 20 minutes of cooling. The heat-cool cycles was ap- 
plied to each specimen six times. The external surfaces of these specimens were then 
examined macroscopically. In addition, about one-third the length of each specimen was 
subsequently mounted longitudinally and polished for metallographic study and for addi- 
tional heat- cool study in a hot-stage microscope. 

The specimens of 0. 40- and of 0. 70- volume-fraction tungsten fiber which were 
studied in a hot-stage microscope were heated from ambient temperature to 1652° F 
(900° C) and back to ambient. The time to heat was about 50 seconds and the time to 
cool was about 10 seconds. The dwell time at temperature was from 5 minutes to as 
long as 2 hours, and the number of heat-cool cycles was from 1 to 18 cycles. These 
specimens were photographed while in the hot-stage microscope. Film speeds of 1 to 
20 frames per second were used. The intent of these particular studies was to deter- 
mine whether damage could be initiated in the specimens by heat-cool cycles and to 
follow the damage in the first cycle, and in some cases in succeeding cycles, by means 
of a photographic record. 

Initial studies which were for the primary purpose of developing methods, were per- 
formed on the transverse (disk) specimens. Subsequent studies were performed on 
longitudinal specimens. The longitudinal specimens were obtained from the 1/3-inch 
(0. 13-cm) lengths which were in turn obtained from the original specimens following 
six heat-cool cycles in the air atmosphere furnace. The 1/3-inch (0. 13-cm) long speci- 
mens which had been mounted, were removed from their mounts and cut into longitudinal 
specimens. The longitudinal specimens were cut to about 0. 25 inch (0. 625 cm) long so 
they would fit into the hot-stage microscope. The transverse and the longitudinal speci- 
mens were examined using conventional microscopy before and after hot-stage micros- 
copy, in order to observe their entire cross sections before and after heating. The area 
under view in the hot-stage microscope was a small portion of the total surface area. 


Results and Discussion 

After the six heat-cool cycles in the air atmosphere furnace, each of the two fiber 
composite specimens sustained macroscopic change or damage. The tungsten fibers, 
which had been flush in both the specimens, were found now to extend beyond the matrix, 
at either end of either specimen. The tungsten fibers also protruded slightly beyond the 
peripheral confines of the matrix in the case of the 0. 70-volume-fraction composite. 

The total length, as well as the diameter, was essentially unchanged in both instances. 
The protrusion was more pronounced at one edge of each specimen, and it was gener- 
ally more pronounced at both edges for the 0. 40- volume-fraction tungsten specimen. 

The edge protrudence of the fibers is shown in figure 22. Further, the 0.40-volume- 








fraction tungsten-fiber composite had an external crack about half-way along its length. 
The crack was partially circumferential and partly longitudinal (fig. 23). 

The protrusion of the tungsten fibers beyond the confines of the matrix may have 
been the result of one or more possible factors, but, because of uncertainties and ques- 
tionable relevancy, this was not pursued. 



Figure 23. - Tungsten-fiber composite after six heat-cool cycles to and from 1600° F 
(877° C) in air furnace. Fiber volume fraction, 0.40. Note external crack. 


The peripheral crack in the matrix of the 0. 40- volume-fraction tungsten-fiber com- 
posite is attributed to a nonuniform distribution of fibers which allowed the generation of 
sufficiently high tensile stresses in the copper on cooling to cause fracture. 

Micro examination of the 0. 70-volume-fraction tungsten specimen indicated micro- 
fracturing of the matrix after the six heat-cool cycles (fig. 24). There was no indication 
of microcracks in the instance of the 0.40-volume-fraction tungsten-fiber- reinforced 
composite. The microfracturing in the 0. 70-volume-fraction tungsten specimen could 
have occurred either while cooling or heating. 

Exposure of a longitudinal section of a 0. 70- volume-fraction tungsten fiber- 
reinforced composite to 1652° F (900° C) for 5 minutes in the hot-stage microscope re- 
sulted in deformation in the matrix, leading to eventual fracturing of the matrix. The 
deformation and subsequent fracturing occurred after the specimen had been at tempera- 
ture for about 2 minutes. The specimen is shown before and after fracture in figure 25. 
The application of an additional 18 cycles resulted in a few more cracks and the accentua- 
tion of existent cracks. Also, matrix recrystallization occurred. This is shown in the 
post-test photomicrograph of figure 26. 
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!a) Strain fracture and recrystallization in matrix. 


(b) Longitudinal cracks in matrix. 

Tungsten-fiber composite after 18 cycles of heating to and cooling from 1652° F 
hot-stage microscope. Fiber volume fraction, 0.70. (Photographed outside hot- 







(b) After heating. Note strain lines in matrix. 

Figure 27. - Tungsten-fiber composite before and after one heat-cool cycle to 1652 F (900 C) 
in hot-stage microscope. Fiber volume fraction, 0.40. (These photomicrographs were made 
from a positive, hence the detail is reversed from previous figure 25.) 
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Although the matrix of the 0. 70- volume-fraction tungsten composite appeared to have 
sustained damage from the previous six heat-cool cycles in the air atmosphere furnace 
and although there was indication of residual matrix porosity (or microdamage from the 
previous six heat-cool cycles), it is evident that the heating cycles in the hot-stage 
microscope produced additional matrix damage. The matrix damage, which occurred 
while at temperature, could have resulted from radial tensile, tangential compressive, 
or shear stresses developed in the matrix upon heating. 

* A 0.40- volume-fraction tungsten specimen was given 14 heat-cool cycles in the hot- 

stage microscope. Matrix deformation was observed on the first cycle during heating 
( (fig. 27). The deformation became more pronounced with each cycle until matrix re- 

crystallization, which occurred at about the 10th cycle, obscured any further effects. 

The calculated elastic stresses that would be generated in a composite with tungsten 
fibers circumscribed by a copper matrix are presented in table VIII. These values are 
approximate for the case where a fair number of fibers are at the surface and do not 
have a circular cross-section. 

In order to assess the significance of the stress values in table VIII, it is necessary 

TABLE Vm. - CALCULATED ELASTIC STRESSES IN TUNGSTEN- FIBER-REINFORCED - 


Reinforce 
ment 
volume 
fraction 

0.40 

0.70 Fiber -139 -956 -139 -956 -139 -956 

Matrix 324 2240 324 2240 324 2240 80 551 

Reinforce 
ment 
volume 
fraction 

0.40 
0.70 



(b) Heating from 80° to 1650° F (26.5° to 1093.5° C) 


Constit- 



Tensile or compressive stress 




c 
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0 

R 

c 

X 

T 0x 


ksi 

MN/m 2 

ksi 

MN/m 2 

ksi 

MN/m 2 

ksi 

MN/m 2 


1 1 3 

78 

1 1 3 

78 

5?, 5 

362 



Matrix 

-48.4 

-334 

11.3 

78 

-42.5 

-294 

71. 1 

490 


79 fi 


79 6 

550 

7Q 6 

550 



Matrix 

-186 

-129 

-186 

-129 

-186 

-129 

66 

455 
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to know the strengths of the matrix and reinforcement of the composite studied. These 
values are as follows: 


Material 

Temperature 

Stress 


°F 

°C 

a 

ult 

b 

T 

(a) 



ksi 

MN/m 2 

ksi 

MN/m 2 

Copper 

80 

26. 5 

19.3 

133 

9.5 

60.9 


1650 

900 

1.5 

10.4 

.75 

5.2 

Tungsten 

80 

26.5 

350 

2410 

— 

— 


1650 

900 

75 

518 

— 

— 


a Copper data based on ref. 20; tungsten data based on 
ref. 10. 


Referring now to table VIII, it is seen that the matrix tensile strengths and, in all prob- 
ability, the matrix compressive and shear strengths can be exceeded in both the 0. 40- 
and the 0. 70- volume fraction composites upon either cooling or heating if there is no 
yielding. The fiber tensile strength could be exceeded in the 0. 70-volume fraction tung- 
sten composite during heating. The stress values in table VIII are based on elastic be- 
havior and might not be reached because of plastic deformation. One would have to make 
a qualitative judgement and say that in the 0.40- volume-fraction tungsten composite, ma- 
trix fracture or deformation could occur either during heating or cooling. Most likely 
however, matrix yielding would occur at least during the initial cycling thus averting 
fracture or unbonding. Fiber damage might possibly occur in the 0. 70-volume-fraction 
tungsten composite during heating, and matrix damage should occur during cooling. 
During heating, the matrix is assumed to be under hydrostatic compressive stress in the 
interior of the composite, but it is not under hydrostatic stress at the surface. This may 
be why the matrix appeared to have fractured while being heated in the hot-stage micro- 
scope; the matrix area viewed was at the surface. The matrix surface area viewed was 
not under hydrostatic stress and was free to experience strain (see fig. 24); the com- 
pressive stresses evidently were great enough to cause fracture. 

The calculated stresses along with the assumption of plastic yielding in low volume- 
fraction composites are consistent with the observed results in that they suggest matrix 
yielding in the 0. 40-volume-fraction tungsten-fiber composite and fiber or matrix frac- 
ture in the 0. 70- volume-fraction tungsten fiber composite. 

The results of the experimental study of the effects of heat-cool cycles on tungsten 
fiber- reinforced copper- matrix composites has shown that merely by heating and cooling, 
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structural changes can be introduced into them. Heating appeared to be more deleterious 
than cooling, but this may have been due to the particular specimen configuration used in 
the hot-stage microscope studies. Thermal gradients may have been introduced into 
these specimens and therefore could have been a contributing factor. But this is not con- 
sidered too likely because of the high thermal conductivity of the copper matrix. One 
can tentatively conclude that high volume-fraction fiber- reinforced composites may sus- 
tain microstructural damage after a few heat-cool cycles or perhaps just as the result of 
cooldown after fabrication. Low volume-fraction fiber- reinforced composites appear 
more able to withstand heat-cool cycles but may ultimately experience thermal fatigue 
because of plastic strains. 


SIGNIFICANCE OF MATHEMATICAL ANALYSIS AND EXPERIMENTAL RESULTS 

Equations for calculating elastic stress and elastic-plastic stresses and strains 
were derived. Calculations of elastic stress were made for a model tungsten-fiber 
80Ni+20Cr matrix composite. Estimates were also made of the maximum elastic -plastic 
stresses and strains that might occur in the model composite. 

The stresses and strains in either constituent was a function of the relative amount 
of that constituent and the calculated elastic stresses could exceed the uniaxial tensile 
strength (and estimated matrix shear strength) of the constituents. By assuming the oc- 
currence of elastic -plastic stresses and strains, rather than elastic behavior, actual 
stresses would be expected to be reduced compared with calculated elastic values and, 
in the case of the matrix, would be reduced to values approximating its yield strength. 
However, the reduction of elastic stresses by plastic strain can result in significant 
amounts of total strain and pose a possible thermal fatigue problem if the composite is 
thermally cycled. 

Experimental studies on 0.40- and 0.70-volume-fraction tungsten fiber- reinforced - 
copper- matrix composites tended to confirm the mathematics. The studies indicated 
that it was possible to generate stresses that were great enough to fracture the matrix 
of the 0. 70- volume-fraction tungsten-fiber composite, presumably upon heating in a 
vacuum to 1652° F (900° C). The stresses in a low volume-percent fiber- reinforced 
composite (i. e. , 0.40-vol. fraction) which was heated to 1652° F (900° C), caused pro- 
nounced deformation of the matrix. Thus, high volume-fraction fiber composites 
(0. 65 < V r < 0. 90) might be prone to fracture of the matrix in a few heat-cool cycles, 
but lower volume-percent fiber reinforced composites may experience matrix thermal 
fatigue after many cycles. However, very low volume-fraction reinforcement compos- 
ites may be prone to reinforcement damage. 
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CONCLUDING REMARKS 


Mathematical analyses, sample calculations, and limited experimentation all have 
indicated that thermal stress, arising from differences in the coefficients of thermal 
expansion of the constituents, can cause structural damage to laminated and fiber- 
reinforced composites when subjected to a thermal transition or to heat- cool cycles. It 
has been postulated throughout that damage is possible as a result of thermally induced 
stresses, and this is felt to have been borne out experimentally; if stress relief can oc- 
cur by plastic flow, preventing immediate structural damage or damage in but a few cy- 
cles, then, under use- conditions involving a multiplicity of heat-cool cycles, thermal 
fatigue could result. The occurrence of thermal fatigue has not been necessarily borne 
out, however, Although this study strongly suggests the existence of a significant prob- 
lem area, there is still very much that has to be understood regarding the susceptibility 
of metallic composite materials to structural and strength damage resulting from ther- 
mally induced stresses or strains due to differing coefficients of thermal expansion. 
Specific composite materials combinations may well possess adequate ductility to pre- 
vent the generation of critically high stresses and simultaneously possess ample thermal 
fatigue strength for heat-cool cycling applications. Additional study, both analytical and 
experimental, appears necessary. 

The problem considered in this study, if of practical significance, might be mini- 
mized by a number of expedients suggested by the mathematics and experimental results; 
these include 

(1) The obvious close matching of thermal expansivities 

(2) Use of materials combinations with good interface bond strengths 

(3) Use of relative amounts of constituents to minimize deleterious effects in the 

crucial constituent, that is, to reduce the stress in the critical constituent by 
increasing its volume fraction 

(4) Use of appropriate spacings (or volume fractions) to reduce shear stress 

(5) Use of intermediate zones of material with intermediate thermal expansivities 

(6) Use of materials with good thermal fatigue properties 

SUMMARY OF RESUL TS AND CONCLUSIONS 

An analytical and experimental study has been conducted to obtain indications of the 
possible effects of differences in coefficients of thermal expansion of the constituents on 
the structural integrity of laminate and fiber- reinforced composites. The study was 
made with reference to metallic systems, specifically to refractory-metal- reinforced 
superalloy- matrix composites. The results and conclusions of this heurestically oriented 
study are as follows: 

1. Approximate equations were derived to permit the calculation of elastic stresses 
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in both laminate and fiber composites containing constituents with different coefficients 
of thermal expansion. The equations were applied to model laminate and fiber compos- 
ites of tungsten- reinforcement and 80 nickel + 20 chromium matrix. Extremely high 
elastic stresses far in excess of the actual or estimated strengths of the constituents 
were calculated for temperature transitions between 80° to 2000° F (26.5° to 1093.5° C). 

2. Approximate elastic-plastic solutions indicated that lesser stresses would occur 
as a result of plastic flow of either constituent. However, relatively substantial amounts 
of strain could be associated with these reduced stresses. 

3. It is postulated that in heat- cool usage, large elastic stresses, sufficient to cause 
fracturing of either constituent or the interface, could quickly occur, on the one hand, 
but, if relieved by plastic deformation, could cause cycle limiting thermal fatigue, on 
the other hand. 

4. Limited experimental studies on laminated tungsten- reinforced - 80 nickel + 

20 chromium matrix composites containing about 50 volume percent of either constituent 
resulted in the following. With a slow heat-cool cycle (to and from 80° and 2000° F 
(26. 5° and 1093. 5° C)), 11 such cycles could cause extensive delamination in a specimen 
with 0. 020-inch (0. 050-cm) tungsten laminae and could cause fracture of the laminae in 
a specimen with 0. 005-inch (0. 0125-cm) tungsten laminae. With a rapid heat-cool cycle 
(to and from 80° and 1600° F (26. 5° and 877° C)) a specimen with 0.020-inch (0.05-cm) 
tungsten laminae could be debonded in one cycle. A specimen with 0. 005-inch (0. 0125- 
cm) tungsten laminae could be severely damaged in three such cycles. And a specimen 
with 0.001-inch (0.0025-cm) tungsten laminae could evidence incipient damage after six 
such cycles. 

6. Limited experimentation on fiber- reinforced composites indicated that micro- 
cracks could be induced in the matrix of a 0. 70- volume-fraction tungsten-fiber - copper- 
matrix composite by a single heating to 1652° F (900° C) in a hot-stage microscope. A 
single heating to 1652° F (900° C) of a 0. 40- volume-fraction tungsten composite in the 
hot-stage microscope produced pronounced matrix deformation. Heating was moderately 
rapid in both instances and may have contributed to the observed effects. Additional heat- 
cool cycles produced some additional fracturing in the 0. 70-volume-fraction tungsten 
composite and a little more pronounced deformation in the 0. 40-volume-fraction tungsten 
composite. Recrystallization of the copper matrix eventually occurred if heat- cool cy- 
cles were continued. Prior to the hot-stage microscope studies, these specimens were 
given six heat-cool cycles to 1600° F (877° C) and the 0. 70- volume-fraction tungsten 
composite suffered matrix microcracking. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 14, 1970, 

129-03. 
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APPENDIX A 


TENSILE AND COMPRESSIVE STRESSES IN LAMINATE COMPOSITES 

The assumptions made for this analysis have already been enumerated on page 5. 
The mechanism by which tensile or compressive stresses are generated have already 
been discussed on page 6. In essence, the resultant strains in either the x or y direc- 
tion (fig. 1) must be compatible and this may be expressed by the following equation: 


“m AT + 


(°x - m y) 
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+ € mxp = “r AT + 


(°X - ro y). 


+ e 


E x E 

m m 


rxp 


(Al) 


But, since a x = a (assumption (4), p. 5), equation (Al) can be expressed as 


%i AT + 


ff xm (1 ' u) 


E 


m 


m a xr (1 ' v) v 
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(A2) 


A second condition that has to be fulfilled is that the total tensile force in the matrix 
be balanced by the total compressive force in the reinforcement and this is expressed by 


P m +P r = ° 


(AS) 


However, the composite consists of N + 1 laminae of matrix material, each of which 
sustains a tensile stress of <x , and N laminae of reinforcing material, each of which 
sustains a compressive stress of a . Equation (A3) can then be expressed as 


( N + Dtcw NWt or A A 

+ r rx = jcn = y a + y CT = 0 

» a A mx A rx m mx r rx 

A c A c c c 


(A4) 


Next, stresses and strains in a complex stress field can be expressed as effective 
or equivalent values by use of the generally accepted Von Mises* relation 
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from which it follows in this instance that 
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°e °x CT y 


(A6) 


The strain in the plastic range may be taken as 
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Since 


°x=S 


it follows that 


e xp e yp 


Assuming incompressibility, it can be shown that 


e zp " 2£ xp 


(A8) 


Consequently, equation (A7) yields the relation 


e xp e yp 2 Ce P 


(A9) 


By rearranging equation (A2) and substituting equations (A4), (A6), and (A9) into it, the 
following equations may be obtained: 
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(All) 


Equations (A10) and (All) can be used to calculate equivalent (or axial) stresses 
when elastic strains only occur in both the matrix and reinforcement or when either one 
or both of the components experience elastic-plastic strain. For the occurrence of 
elastic strains only, the plastic strain terms in equation (A10) and (All) are omitted. 
Equivalent stress and strain have been used throughout to obtain indications of maximum 
strains . 

Elastic-plastic solution. - For elastic -plastic solutions, however, equations (A10) 
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and (All) must be solved by trial and error or by iteration using the stress- strain curves 
of the materials. This is relatively rapid when using computers. However, for the pur- 
pose of obtaining magnitudes, an expedient graphical method, although not necessarily 
rigorous, has been used in this report. The graphical method is based on the simulta- 
neous satisfaction of two requirements: (1) the strains in the two constituents must be 
compatible and (2) the forces must be in equilibrium. 

First, for compatibility of strain in the x (or y) direction the following condition 
must be fulfilled 


e tmx + a m 


AT 


€ trx + “r AT 


(A12) 


However, for simplification it is assumed that all the strain is plastic. This assumption 
allows the use of Von Mises' equation for effective plastic strain, that is, 
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ep 
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Further, using the usual assumption of incompressibility, it can be determined that 


e ep = 2e xp = 2e yp 


(A13) 


Thus, using equations (A12) and (A13) and the assumption of complete plastic strain 


e - e « 2(a - o' )AT 
emp erp v r nr 


(A14) 


Next, forces must be in equilibrium and satisfy the following equation 


V a + V (7 =0 

m em r er 


(A15) 


Equations (A14) and (A15) can be solved using the appropriate stress-strain dia- 
grams for the matrix and for the reinforcement. Corresponding values of o em , e epm , 
a gr , and e epr can be obtained from these diagrams, which simultaneously satisfy equa- 
tions (A14) and (A15). It should be noted that the values of a and e on any stress- 

r 

strain diagram are assumed to correspond to values of a g and e gp . The values thus 
obtained will be estimates of the elastic-plastic stresses and strains and will be more 
accurate when the elastic strain is small compared with the plastic strain. 

Another point to mention is that stress -strain diagrams for the temperature ex- 
tremes (i. e. , 80° F (26. 5° C) or 2000° F (1093.5° C)) have been used, which assumes 
the absence of any stress relaxation. 
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APPENDIX B 


ELASTIC SHEAR STRESSES IN LAMINATE COMPOSITES 

The assumptions made have already been enumerated on page 7, as was the mecha- 
nism of stress generation. The region of interest is shown in figure 1, and the element 
of interest is located at y = y Q and is near the tip of the composite (i. e. , x « 0). A 
typical element is shown in figure 3(a). 

The governing equation is derived from equilibrium of forces in the x direction. 
From figure 3(b), 
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From which it follows 


(Bl) 


da 


rx 


dx 


- 2t = 0 
mzx 


(B2) 


The matrix is under pure shear as illustrated in figure 4. From the geometry of fig- 
ure 4, the definition of shear strain, and Hooke's law, the following relations can be 
written, letting e x equal the displacement at the tip of the composite 
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from which it follows 
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Also from Hooke's law 
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From the assumption that 


it follows that 


a = cr 

x y 


e y = e x 


which, when substituted into (B6a), results in 
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But since 


e = 
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du 

dx 


it follows that, for the reinforcement, 
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Using equations (Bl), (B5), and (B8) and dividing by t results in 
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or 


d 2 u 4G m^ 1 “ ^ r u 4G ^ ' ^r e av x 


dx 
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Since equation (B9) is linear, it can be solved by superposition, 
neous solution, 

u h =Ae- X l X/t ” t Be A > X/t -" 

where 
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The particular solution yields 


and the solution now becomes 
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where is given by equation (Bll). 

The boundary conditions for equation (B12) are as x — °o 


d < u - e av x > 
dx 


- 0 


and at x = 0 
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Using these conditions yields 


(B9) 

For the homoge- 

(B10) 

(Bll) 
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. *l x Am A.,/t X«xA 
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(B14) 


Applying the second condition 
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and u now becomes 
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The shear stress is obtained by substitution of equation (B14) in equation (B5) yielding 
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The average strain is obtained from the condition that at x » t (away from a tip), 
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(B17a) 
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It is also apparent that average stresses may be equated giving 
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Equating (B17a) and (B17b) and using equation (B18) yield the relation 
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and the average strain e ay , becomes 
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Using equation (B20) in equation (B16) yields the equation for shear stress at the inter- 
face or in the matrix at any point x 
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The limitations of equation (B21) are discussed in appendix F. 


(B22) 
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APPENDIX C 


TENSILE AND COMPRESSIVE STRESSES IN HIGH-VOLUME-FRACTION, 
FULL-LENGTH FIBER-REINFORCED COMPOSITES 


The assumptions to be used in this derivation have already been enumerated on 
page 29. The mechanism by which stresses are generated have already been alluded to 
on page 29 and is essentially as follows. After cooling, the average strains at a distance 
away from the ends must be compatible (i. e. , the total strain in the fiber must be equal 
to the total strain in the matrix. And the forces must be in equilibrium. There will be 
equilibrium radial, tangential, and longitudinal stresses in both matrix and reinforce- 
ment (i. e. , Op^, cjg, and ct x ). First, consider the longitudinal or x direction stresses. 
Since the strains in this direction must be the same for both fiber and matrix; that is, 
the change in length of the fiber must be equal to the change in length of the matrix, the 
following relation between strains may be written: 
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or from assumption (3) (p. 29) 


(C2) 


Since there must be equilibrium of forces, the following condition must also be satisfied 


P + P = 0 
m r 


(C3) 


Each fiber of area is assumed to be interacting with an area of matrix a m , and 
since there are N fibers, equation (C3) can be represented as 


Na cr + Na cr =0 
r r mm 


(C4) 


Dividing through by the area of the composite A c , gives 


Na r a r ' Na m°m _ Q 


(C5) 
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However, since area fraction and volume fraction are equal, equation (C5) can be written 
as 


V a + V a = 0 
r r mm 


(C6) 


Since there can be no plastic deformation in a condition of hydrostatic stress, the terms 
for plastic deformation are zero in equation (C2). Equation (C2) can now be written as 


a m AT + (1 - 2v) = AT + (1 - 2v) 
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m 


E 
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Equation (C6) can be substituted into equation (C7), and by algebraic rearrangement the 
following equation representing the longitudinal stress in the matrix of the composite can 
be obtained 
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In a similar fashion the following equation representing the longitudinal stress in the 
reinforcement (fiber) can be obtained 


V m E r E m( a r ' “m> iT 
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(C9) 


Since a = = a D has been assumed, equation (C8) can also be used to obtain the 

X u rt 

radial and tangential stresses in the matrix and equation (C9) can also be used to obtain 
the radial and tangential stresses in the reinforcing fiber. Obviously, equivalent 
stresses and strains will be zero, but not necessarily the principal stresses. 
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APPENDIX D 


TENSILE AND COMPRESSIVE STRESSES IN LOW- VOLUME-FRACTION, 
FULL-LENGTH FIBER-REINFORCED COMPOSITES 


The assumptions made for this analysis have already been enumerated on page 30 
and, the mechanism whereby the stresses are generated has been discussed on page 30. 
The basic equations that will be used have been derived by Timoshenko (ref. 21), and are 
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where a, b, and c are radii of the cylinders shown in figure 13(d), p^ and p Q are in- 
ternal and external pressures, respectively, and 6 is the difference between fiber di- 
ameter and matrix element inner diameter after heating or cooling (in the absence of any 
mutual restraint). 

Equations (Dl) and (D2) may be used to express the stress in the matrix in which 
case p Q is assumed to be zero and the inner and outer radii are taken as b and c, re- 
spectively. The resulting equations are 
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Now, if b and r, are divided by c to make the equations dimensionless, the following 
equations result: 


where 
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Further, if 


r is made equal to b, the inner radius, the preceding equations then become 
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and give the maximum stress which is at the fiber interface. Equations (Cl) and (C2) 
may also be used to express the stresses in the fiber. In the case of the fiber, the inner 
radius a becomes zero and the internal pressure is zero. However, the fiber is under 
some external pressure p Q . In a manner analogous to that just used, dimensionless 
equations may be derived for the stresses in the fiber and these equations are 

a rR = -P 0 < D10 > 

"re = -Po < Dll > 


Equations (D8) and (D9) give the stress at the interface between the fiber and matrix, and 
equations (DIO) and (Dll) indicate constant stresses anywhere along the radius of the 
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reinforcing fiber. Heating will reverse the stress directions. 

It is apparent that the external pressure p Q on the fiber must equal this internal 
pressure on the matrix; that is, 

P Q = p i = p ( D12 ) 

The value of p can be obtained from equation (C3) by setting 6 equal to the difference 
between the fiber outer radius and the matrix inner radius after heating or cooling; that 
is, 


5 = b ( Qf m " a r ) AT 

Further in the derivation of the equation for p, the fiber and matrix are assumed to have 
different elastic moduli but equal values of Poisson's ratio, and the equation for p may 
be made dimensionless by dividing c by b. The resultant equation is then 


<“m ' a r> iT E m E r< 1 ‘ ^ 

E r [ 1 e 2 a - V) + l + I/)] + E m (i - (?)<i - v) 


(D13) 


Equations (D8) to (Dll) and (D13) can be used to calculate the elastic stress in the 
matrix and in the fibers in principal directions transverse to the length of the compos- 
ite. The preceding equations will give conservative results since they are predicated on 
the assumption of no stress interaction between adjacent matrix cylinders. Actually, 
there are stress interactions, and the strain at the periphery of each matrix cylinder 
would approach zero because of them. Reevaluation of the constants in equations (Dl) to 
(D3) with the assumption that strain equals zero at c might be a better assumption to 
make particularly when there is a considerable amount of fiber; hence, equations (Dl)to 
(D3) are probably conservative. 

The longitudinal stresses in the reinforcing fiber and in the matrix may be derived 
by assuming compatibility of strain and equilibrium of stress. Doing this results in the 
following equations: 
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Inspection of equations (D6) and (D7) indicate that for the matrix 




(D16) 


Inspection of equations (DIO) and (Dll) indicate that for the fiber 


( a e + °r)„ = - 2 P 


(D17) 


By substitution of equations (D15) to (D17) into equation (D14) and by algebraic rearrange- 
ment, the following equations may be derived: 
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Equations (D18) and (D19) may be used to estimate the elastic stresses in the fiber 
and in the matrix. 


Elastic-Plastic Stresses and Strains 

Equations (D8) to (Dll) and (D19) are predicated on the assumption of elastic defor- 
mation and contain no terms for plastic deformation. But they can be used, in conjunc- 
tion with the strain invariance principle, to estimate elastic-plastic stresses and strains. 
The strain-invariance principle (ref. 12) asserts that, in thermal strain problems, the 
total strain, that is, the elastic plus plastic strain, can be estimated by assuming com- 
pletely elastic behavior and then calculating the elastic strain associated with a calcu- 
lated elastic stress, even though this stress be beyond the yield strength. 

The previously mentioned equations will be used to calculate elastic stress values 
for the 6, R, and x directions. Using these calculated values along with the strain in- 
variance principle, the total strain in each of these directions will be calculated for 
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fiber- reinforced composites with fiber fractions up to 0.65. Equivalent total strains will 
be calculated for a given set of principal strains using the relation for e ^ in equation 
(35) of reference 22, that is, 


e 



e 





1/2 


Having calculated e ^ for a given value of V r (or V m ), the appropriate stress-strain 
diagram (fig. 7) will be used to obtain estimates of the corresponding values of e and 

©P 

Og. With the value of known, the following equation (i. e. , eq. (5) of ref. 22) 
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will be used along with the stress-strain diagrams to get corresponding values of e 
e , and a . The values of e and a correspond to e and a of the appropriate 

c cp c p 

stress-strain diagram. Values of e e ^., e g p, and o e , which satisfy both equation (D20) 
and the appropriate stress-strain diagram, are first approximations of equivalent 
stresses and strains and will indicate the plastic strain as well as the total strain that 
has occurred and what the corresponding stress is. 

Stress-strain diagrams for the limiting temperatures (80° F (265° C) and 2000° F 
(1093. 5° C)) will be used. The diagrams are used on the assumption that no relaxation 
occurs during the temperature change. 
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APPENDIX E 


ELASTIC SHEAR STRESS IN FULL-LENGTH FIBER-REINFORCED COMPOSITES 


The assumptions to be used in this analysis as well as the mechanism which gener- 
ates shear stress in the matrix have been discussed on page 27. 

The strain in the longitudinal or x direction of the composite may be expressed as 
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Referring now to equations (D13) and (D17), it may be seen that, for a given composite 
having a specific relative amount of reinforcement and exposed to a specific AT, 
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Equation (D19) can also be evaluated for the same set of specific conditions as equation 
(E2), so that for a set of specific values 
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By placing equation (E3) into equation (El), the following equation may be obtained 
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Figures 14 and 15 show that equilibrium of forces between the fiber and its associ- 
ated cylinder of matrix may be written as 
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dx - 
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The equation for shear stress in the matrix of the fiber-reinforced composite is deter- 
mined, in an analogous manner to that used in appendix B, for the laminate matrix. 
Figure 4 can be used again except 1/2 t is now (D m - D r )/2 and the shear stress is 
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Substitution of (E8) in (E7) results in 
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Next, differentiating equation (E5) and substituting it into equation (E7) result in 
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An alternate of equation (Ell) is 


d 2 u 8G m (l-,K) rl , 8G m (l-K, r )€ avX 


,2 (D -D )D E 
dx v m r ; r r 


(D - D )D E 
'm r ; rr 


(E12) 


Solving this equation as before results in 
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(E13) 


-A 2 x/( D m - D r) 

u = e ay x + Ae + Be 


where 


2 8G m^ - K ^ D m - D r> 


D, 


Applying the same limit conditions as previously allows evaluation of the 
and results in 

< D m - D r> £ av „V/< D m - D r> 

u= £ a v x + : e 


Equation (E15) is substituted into equation (E8) resulting in 

-X 2 x/(D m - D r) 


2G m< D m - D r> e av e 


m0x 


< D m - D r> x 2 


or 


2G e e 
mav 


-V/< D m- D r> 


m0x 


Since strains must be compatible, it follows that 


- = (1 - vK r ) ^ + a r AT = (l - ^ + « m AT 
E„ 


-av,x '-r' R ~r " m E ~m 


m 


It also follows that 


From which it follows that 


D k + ( D m ~ D IK = 0 


d; 


a. 


mx 


D 2 - D 2 
m r 


(E14) 

constants 

(E15) 


(E16) 


(E17) 


(E18) 


(E19) 


(E20) 
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Using the relations of equations (El 8) and (E19), it follows that 


rx 


(i - lag (i - Ku m ) d; 


E_ 


E„ 


m 


b 2 -D 2 
m r 


= (« m - « r )AT 


(E21) 


This equation may be arranged to give 


(a - ffjAT 
' m r' 


rx 


(1 - Kv r ) (1 - K v ) D* 
+ 


E. 


E m (D 2 - D 2 
m r 


(E22) 


or expressed in terms of strain 


(a -a )AT 
' m r' 


'av, x 


( 1 - K ^/ E r\ D r 


1 + 


(1 - K ^r) \ m/ - D 2 


(E23) 


This equation may be substituted into equation (El 7) to give the final equation which ex- 
presses the shear stress in the matrix of the fiber-reinforced composite 


where 


7 m0x 


(a 
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- ff r )AT 2G m e 


-A 2 x/( D m -D r ) 


1 + 


< 1 - 

(1 - Ku r ) 




8G (1 - vK)(D m - D ) 
m 'r m r 7 


D 


(E24) 


(E25) 


The relative breadth of matrix is subject to the limitations discussed in appendix F. 
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APPENDIX F 


ASSUMPTION INHERENT IN SHEAR-STRESS ANALYSES 

It is assumed that uniform shear will take place in the matrix as long as the matrix 
thickness does not exceed the reinforcement thickness in the case of laminate compos- 
ites or as long as the interfiber spacing (i. e. , matrix thickness) does not exceed the 
fiber diameter in the case of fiber reinforced composites. 

The shear lag analysis then applies to laminate composites whose reinforcement 
volume fraction is greater than 0. 5; that is, V r > 0. 5. 

In the case of a fiber- reinforced composite with a hexagonal array of fibers (see 
sketch) 



it can be determined that the ratio of fiber to matrix is 


V r = 


1 nd l 
3x - —L 

6 4 




Equation (Fl) can be reduced to 


7Td* 


V r = 


2 X i(2 df ) 2 ii Vi 


0.23 


when 5 = d^ 
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(Fl) 


(F2) 
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Therefore, the shear lag analysis, under the assumptions made, is applicable to 
fiber-reinforced composites, with hexagonal arrays of fibers and with 

V r >0.23 

It can be readily shown that for square arrays of fibers 

V r >0. 19 

As a generalization, the shear lag analysis is assumed valid for fiber- reinforced 
composites when the volume fraction of fiber is in excess of about 20 percent. With the 
exception perhaps of whisker-reinforced composites, most filament- reinforced compos- 
ites will have a V r greater than about 0.25 to 0.30; hence, the shear analysis is appli- 
cable for most practical fiber composites. In the case of laminates, it is conceivable 
that V r < 0. 5, in which case equations (B3) and (B2) are not necessarily applicable for 
the lower amounts of reinforcement. 

The limitation of the shear lag analysis is assumed to be that the matrix thickness 
does not exceed the reinforcement thickness as illustrated in figure 28(a). If the matrix 
thickness exceeds the reinforcement thickness as shown in figure 28(b), it is then as- 
sumed that the shear deformation occurs in a region not exceeding r/2 in thickness. 

On this assumption the actual shear angle at the tip y g (of fig. 28(b)) would be greater 
than the shear angle assumed by use of the shear lag analysis y g. Hence, the tip shear 
stresses associated with y^ would be greater than those associated with y^, and shear 
stresses based on a shear lag analysis would be lower bound values where m > r. The 
derived equations (B21) and (E24) will give lower boundary shear stress values for those 
cases where V r <t 0.5 and V r < 0.2 for laminate and fiber composites, respectively. 
A more complex analysis such as that presented in reference 23, would be necessary to 
treat the problem of lower volume-fraction reinforced composites. 
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Matrix 
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Reinforcement^ 



(a) Assumption made in 
shear lag analysis; 
that is, matrix spacing 
is equal to or less than 
reinforcement thickness 
or diameter and shear 
in the matrix extends 
for a distance r/2. 
Therefore, yj = 2u/r. 


(b) Matrix spacing greater than 
reinforcement thickness or di- 
ameter. When shear in the 
matrix occurs for a distance of 
r/2 rather than m/2, y 2 ■ 
(2u/m) <y 3 ■ 2u/r. 


Figure 28. - Shear occuring in matrix whose thickness is (a) equal 
to or less than the thickness or diameter of the reinforcement or 
(b) greater than the thickness or diameter of the reinforcement 
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